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We demonstrate the relation of the infrared anomaly of conformal field theory with entropy 
\N ■ considerations of finite temperature thermodynamics for the 3-state Potts chain. We compute 

the free energy and compute the low temperature specific heat for both the ferromagnetic and 



X 



Abstract 



> 

^b ' anti-ferromagnetic spin chains, and find the central charges for both. 

1 Introduction 
(N 

The conformal field theory treatment of quantum spin chains at positive temperature T deals 

with systems of size M in the limit 

^- M^oo, T^O, TM fixed. (1.1) 



This limit is discussed in terms of a variable q = exp(2iTv/MT) where v is the speed of sound. 
The modular invariant partition function is computed in terms of q, and one of the important 
results [||, ^ is that as g — > 1, the free energy per site is given as 

f = eGs-^T' + o{T') (1.2) 

ov 

where c is the central charge as determined from finite size corrections to the ground state energy 

EGs = Meo- — + o{M^). (1.3) 

However the limit ( |l.lD is not the limit of thermodynamics. This limit is 

M ^ oo, T fixed. (1.4) 
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Here the low temperature behavior of the specific heat is obtained by letting T ^ after the 
limit ( |1.4| ) is taken. 

These two limits will give the same result if there are no additional length scales in the problem. 
In this case the result is obtained P, Q that the specific heat C is 

C^—T. (1.5) 

For the anti-ferromagnetic 3-state Potts chain, the low lying order one excitations in the 
limit (1J-), were used in |^ to compute the partition function, and the result ( |l.5D was obtained. 



This was accomplished by adding up the order one excitations found from the Bethe's equations 
for the model [Q, |5|. The partition function obtained in this manner is the modular invariant 
partition function of conformal field theory. 

The counting of states in [Q] depends on the fact that the momenta of order one excitations 
obey not only a fermi exclusion rule Pj 7^ Pk, but also have additional exclusion rules: The number 
of states near e(P) = diminishes as the number of order one excitations in the system increases. 
Correct counting of the states, incorporating these exclusion rules, gives a partition function 
which has a central charge smaller than that of fermions. In the case of the anti-ferromagnetic 
3-state Potts model, there are three kinds of quasi-particle excitations, which, for purely fermionic 
exclusion rules, would give a central charge of 3/2, whereas the actual central charge of the model 
is 1. 

In this paper, we use the thermodynamic limit ( |1.4D to obtain the low temperature specific 
heat of the 3-state Potts chain from Bethe's equations, using the methods of pi 1^, ^. In section 
2, we write the Bethe equations for the finite lattice and introduce the completeness rules [^ for 
(5 = 0. We find it convenient to study the ferromagnetic and antiferromagnetic cases using two 
different sets of integral equations. In section 3, we write the free energy for the ferromagnetic 
case in terms of only one integral equation, and compute the linear term in the low temperature 
specific heat. In section 4, we do the same for the anti-ferromagnetic chain, this time in terms of 
two integral equations. We obtain the central charge of the conformal limit of both spin chains 
from the linear term in the specific heat, using equation ( |1.5| ). In section 5, we consider the sector 
Q = l. 

In section 6, we discuss how the counting of states in the finite size system [^ is incorporated 
into the discussion in terms of densities in the thermodynamic limit. We find that in the ther- 
modynamic limit, the counting of states is contained in the entropy, written in terms of densities. 
The densities are related through the thermodynamic limit of the Bethe equations ( |3.1[ ), ([4.1|), 
and this relationship corresponds to the way the number of available states depends on the num- 
ber of excitations in the finite size system. These equations become particularly important when 
evaluating the low temperature specific heat, where we need to consider the equations in precisely 
the limit which corresponds to P ^ in order to extract the linear term in the specific heat. 



It is important to note that the 3-state Potts model is the D4 model in the classification of 
Pasquier |1C], which, by orbifold construction [^, 12 1, has certain sectors of eigenvalues which 



overlap with the A^ RSOS model. From this construction, it is to be expected that the ther- 
modynamic quantities of the two models are equal as long as the sectors which dominate the 
thermodynamics are common to the sectors which overlap. Indeed, the ground states of both 
the ferromagnetic and anti-ferromagnetic ends of the D4 model are the same as the critical A^ 
model at the boundaries of the III / IV and //// regimes respectively. Thus the central charges 



computed for A^ |13, 14, 15 1 coincide with the central charges of the D4 model [1£, 17]. We further 



note that the classical two dimensional anti-ferromagnetic 3-state Potts model which is critical at 



T = is equivalent ||l^, |l^, Q to the 3-coloring problem and also has central charge c = 1 1 24 ] . 
The thermodynamics of the A^ series was studied in ||l5|] , where the integral equations for the 
free energy, the central charges and order one excitations above the ground state were found. The 
thermodynamics of other affine Lie algebras were also studied in ||2^. In this paper, however, we 
obtain different sets of integral equations, which display a more direct relationship to the modular 
invariant partition function discussed in [^ for the anti-ferromagnetic chain and in [^3| for the 
ferromagnetic chain. The fact that there are different sets of integral equations for the model is 
related to the fact that the modular invariant partition function can be expressed as sums over 
different sets of quasi-particle excitations |2^ 

2 Formulation 

The three state Potts hamiltonian is 

2 '-' 



H = ±^ Ei^. + ^j + Zj^U + 4^^-+i} (2-1) 



with periodic boundary conditions, where M is the number of sites in the chain, and the matrices 



Xj and Zj are 

Xj = I (g) I ... Xjth (g) . . . (g) /, 

Zj = I®I® ...®Zjth® ...®I, (2.2) 

/ being the 3x3 identity matrix and X and Z are 3x3 matrices with entries: 

Xij = (5ij+i(mod 3), Zij = Si^ju'-^ , uj = e^^'/^ . (2.3) 

The hamiltonian with the (-I-) — sign is referred to as the (anti-)ferromagnetic Potts chain It 
commutes with the spin rotation operator, whose eigenvalue is e^'^*'^'^, Q = 0, ±1. The eigenvalues 



of the hamiltonian ( |2.1| ) are derived from functional equations |13, 15, ^, 26, 27| from which we 
find §: 

IT / l\/l 

(2.4) 



E = J2(^ot{iXj + ^)-^, L = 2{M-\Q\), g = 0,1,2 



where the set {Xj} satisfy the Bethe equations: 



■sinh(i7r/12 — Xj 



sinh(i7r/12 + Aj 



2M 



(_l)M+i TT sinh(i7r/3 - {Xj - Xk)) 
^ ' ^}^smh{i7r/3 + {X,-Xk)y 



1,...,L 



(2.5) 



Not all solutions of the equations (^^) correspond to eigenvalues of the hamiltonian (|2.1| ) . The 
equations do not impose sufficient restrictions on the set {Xj}. There are additional conditions, 
which ensure that the energy is real, postulated from finite size studies in Q, where the spectrum 
of H was classified. We introduce these conditions by writing ( p.5| ) in logarithmic form. The 
solutions of (2.5) which correspond to eigenstates of the hamiltonian ( p.l| ) fall into 5 classes. 



where below Xj is a real number: 



A 



Xj, Xj 
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A, + -, 



V 



\2s 



^7^ 



A7''' 



^■-T- 



V 



\ns 



ITT 

A. ± -. 



V 



(2.6) 



The last three always occur in complex conjugate pairs. Since we are interested in the infinite 
lattice limit, the imaginary parts are assumed to be exact. We denote the number of each type 
of root a £ {+, — , 2s, —2s, ns} by rua- 



We rewrite the Bethe equations (2^) to explicitly display the different types of roots (^ 



Let 



Then (p.5|) become, for a = + or — : 



hiX) 



sinh(f - A) 
sinh(^ + A)' 



(2.7) 



-1) 



M+l 



sinh( 



§-Apn2M 



sinh(^ + A^) 



m^ 



m2s 



n MA° - A+) n MA," - AT) n kxj - xfnx^ - xf*) 



'J 



k=l 



fc=l 



fc=l 



n MA^ - A-^^)/.(A^ - A-^-) n MAJ - Ar)/.(AJ - A5 



ns*\ 



(2.^ 



k=l 



k=l 



whereas for /3 = 2s, —2s or ns, the equations for each complex conjugate pair are multiplied 
together: 



sinh(f| - A^") sinh(f| - Af )n 2M 



sinh(§ + A°) sinh(§ + Af ) 



rri-^ 



n h{x^ - xtMxr - A+) n ma,- - A,-)MAr - a,-) 



fc=i 



k=l 



n h{x^ - \f)h{\f - \f)h{\'^ - xinnxr - >^in 



k=l 

m-2s 



n h{\^ - \f')h{\f - \f')h{\^ - xfnKxr - K 



■2s*\ 



k=l 

rrlna 



n KX'^ - xr)h{xf - xr)h{x^ - XDKX^* - XD 



(2.9) 



fc=i 



We follow ref. Q in taking the logarithm of equations ( |2.8[ ), ( |2.9D . To do this, we define the 
functions ta and Qa,i3 in the following way: 



where f±2s = -I fns 



-2iln 



sinh{i7r/12-A=^) 



-2i/aln 



"sinh(J7r/12+A*) 
sinh(i7r/12-A") sinh(i7r/12-A"*) 



sinh(J7r/12+A") sinh(i7r/12+A°'*) 



for a = lb; 

for a = ±2s,ns 



(2.10) 



-1/2 and f± = 1. For a,/3 = ±: 



i^\ 



-iln 



a/3/i(A^" - Af ) 



e«MA," - A^ = - 

For a = =b, /3 = ±2s, ns or a = ib2s, ns /3 = it, 

0a/3(A," - Af ) = -i/„ ln[e„^/i(A^" - Af )/i(A- - Af ) 
where e+,_2s = e-,2s = —1, = 1 otherwise. For a,/3 = ±2s,ns 

e«^(A," - Af) = -iMnL^/.(A," - X^MX'} - Af )MAr - Af)/.(Af - X^) 



(2.111 



(2.12) 



(2.13) 



where e2s,2s = e-2s,-2s = — !> 1 otherwise. (Note that in [^, the functions i^ and 6a/3 for 
a = ziz2s,ns were defined without the factor /q,. This will change the completeness rules some- 
what from those presented in y], but is necessary in order to have positive densities.) Here, all 
logarithms are taken so that |Im lnz\ < vr, and the functions ta and 0^^ are defined so that 



ta{Xj J 



if Re (A") = 0, 



e„/3(A°-Af) = if Re(A°)=Re(Af). 
The logarithmic Bethe equations are written in terms of these functions: 



1 



Z(Xi) = ^ = —tJX'' 



1 



m/3 



M 27r 



27rM 



E E0"/3(\"-^f)' 



(2.14) 



(2.15) 



P=±,±2s,nsk=l 



where I'j are (half-) integers. We now present the completeness rules for the integers /^. It is 
only necessary at this point to discuss the completeness rules for (5 = 0. It will be shown in 
section 5 that the Q = ±1 sectors are identical to this sector in the thermodynamic limit. The 
completeness rules of H for Q = in the notation introduced here become: 

1. /^ and /^** are distinct (half-)integers, are chosen from the same set of m^ + m2s integers, 
and I^ = /|'* ^ , where h represent a "hole" or missing integer. Therefore the set {lt^+{l'j^^ 
fills the interval — l/2(m+ + m,2s) to 1/2(?tt-+ + m2s)- 

2. I^ and /^ * are distinct (half-)integers, are chosen from the same set of ?n,_ + ?n,_2s integers, 
and /^ = 11"^- Again, the set {I~} + {Ij'^''} fills the interval. 

3. /^'^ are distinct (half-)integers chosen from a set of 2m_ + 2m_2s + f^ns (half-)integers. 

4. The spacing between "available" integers, the set of integers {/^ } + {I^ ^}, is 1. 

We see that + integers correspond to missing 2s integers, and the same for — and —2s. In addition 
to these rules, there is a sum rule for tti^: 

m+ = 2nns + 3m_ + 4m_2s 
mns + 2m„s + 3m_2s + 2m_ = M (2.16) 

This sum rule is responsible for restricting the maximum integers I^ax ^s a function of the number 
of excitations in the system. 

We make the assumption at large M that the rules (1) and (2) imply the equality 

A+ = Af ^ AT = Xf ^ (2.17) 

This appears to be true from numerical results, has been proven for order one excitations |5|, and 
shown to be consistent for all excitation densities in the thermodynamic limit. 

We now take the thermodynamic limit M — > oo of the Bethe equations ( |2.15| ), with A fixed. 
When we do this, we lose the information contained in the rules (l)-(3) about the maximum 



integers. We rewrite the functions Qa,p and ta in terms of the real part A^, using (2.6), and take 
the derivative of Z[\) with respect to A in the thermodynamic limit. We obtain the following set 
of equations: 



.±<^^ - ^r.± - 1 r 



PfW = -K/i2W-^[K^/sHp;-p;n+K^/s*iP;-Pp) + {K/i2 + K/i2}*P: 



pf'w = hK^/uW-K^/,ix))-^[K^/,*{p;-P,'n+K^/s*iPp-pp 



+ {K/12+K/12}*P: 



prw 



where 



27r '"/^ '"/^ 27r 

+ (^W3 + K 



2 (-^7r/12 



^7r/3 "T ^^7t/3' * Pp 



Pt 



lim 



A/-.00 M(A/^.+i 

the convolution * is defined as: 



A/, 



+ ^./12)*(/'p 



Mm 



Pv ^' + Pv 



Pf) 



/oo 
-oo 



and the kernels K^{X) are 



KtW 



±2 sin 2a 



(2.18) 
(2.19) 

(2.20) 
(2.21) 



cosh 2A =F cos 2a 

In writing equations ( 2.1§| ), we did not make use of the relationship between holes and parti- 
cles ( |2.17 ), which imposes a relationship between the densities in equations ( ^.18 ). The assump- 
tion ( 2.17| ) implies that when particle integers of, say, + are equal to the "hole" integers of 2s, 
their corresponding rapidities are are equal. Therefore, in light of the definitions ( p. 19 ), the total 
densities of + and 2s are equal (and those of — and —2s as well), and the particle densities are 
related in a simple way: 



-2s 



±2s 



p+(A) = pf(A), p,-(A)=pr'^(A), pt'%X)=pt{X)-pt{X) 



(2.22) 



This allows us to rewrite the density equations ( |2.18 ), in terms of three independent particle 
densities. It is convenient for further computation do this separately for the ferromagnetic and 
anti-ferromagnetic spin chains. 

3 Ferromagnetic chain 

For the ferromagnetic hamiltonian, we know from |^, ^ that the order one excitations are -|-, — 
and ns. We therefore choose to rewrite (p. 18) using (2.22), as 



ptW 

PT{X) 

pTW 



6 



+ Ki*(p+ + p;)-i^2*pr 



vr cosh 6A 

Ki*{p^ + p-)-K2*p';^ 

K2*{Pp +Pp), 



where the kernels are 



Ki{X) 



18 A 
TT^ sinh 6A ' 



i^2(A) 



vr cosh 6A 



(3.1) 
(3.2) 



The particle densities in equations ( |3.1| ) above are now all independent of each other, there are 
no additional constraints. 

In the thermodynamic limit, the sum rule ( |2.16| ) becomes a relationship between total particle 
densities £)„, 



rrin 



D^ = hm -^ = / dXp^. (3.3) 

However, we find we do not need to impose the sum rule as an additional restriction on the 
densities in ( |3.1| ), as it is contained in those equations already. To see this, we take the fourier 
transform of the first two equations in (^]^) and evaluate at k = 0. This gives exactly the 
relationship ( |2.16| ) divided by M. 

In Q, the sum rules (2.16) were found to give rise to the infrared anomaly, that is, to the 
diminishing of the number of states near P = for the anti-ferromagnetic case, and thus to 
exclusion rules beyond those of fermions. Here, although we lose information about how the 
maximum integers change as a function of m-Q, when we take the thermodynamic limit, we still 
retain a restriction between the densities which contains some of this information. This restriction 
will allow us to retain the concept of correct counting of states in the thermodynamic limit. 

The free energy is 

F = E-TS (3.4) 

evaluated at the stationary point with respect to independent particle densities, where S is the 
entropy of a state with fixed densities p^ and E is the total energy of the state, for large M, the 
entropy is: 

S = M Y. r dX(pflnpf-ppnp^-pfilnpt), (3.5) 



where ph = pt — Pp- The energy E is the thermodynamic limit of equation ( |2.4| 

2M 

where e"(A) is the energy associated with a root of type (5: 



r 9/vf 

E= Y. J d\p'^{X)e^{X)--^ (3.6) 



+ ,,, ±l-2isinh2A „,,,, -2^/3 - 4isinh4A 

e (A) = ^, e (A) = — , 

^ ^ 2cosh2A^^/3 ^ ^ l + 2cosh4A 

e^'^^'iX) = ^l~^^sinh2A ifcl-isinh2A 
2cosh2A=F\/3 cosh2A 



The energy in ( |3.6D is not manifestly real. However, using (]^^, we find that the energy can be 
re-expressed in terms of only the independent particle densities, and depends only on p^: 

^=-7f + /-co''jSih6A''?W- (3.8) 



This expression for the energy is manifestly real. Minimizing 
particle densities p^,p~,Pp^, we obtain the free energy per site, 

6 



f = eo-T 



dX , , 
oo TT cosh 6A 



|) with respect to the three 
ln(l + e-^''W/^), (3.9) 



where 



eo 



2 

7! 



°° 6e2«(A) 
-oo 7rcosh6A 



4 
vr 



3\/3 



-2.81284. 



(3.10) 



and the densities e^ are defined as pf^/Pp = exp (e^/T), and satisfy the nonlinear integral equa- 
tions: 

6 



6+(A) 
6- (A) 

e""(A) 



cosh 6A 

e+(A) 



Ki * In 



(1 + e-^-'/^)(l + e-^"/^) + i^2 * ln(l + 6"^"^^) 



6 



(3.11) 



cosh 6A 
TK2 * [ln(l + e-^^('^)/^)(l + e~^"(^)/'^) 

The functions e" are also referred to as dressed energies. Note that (3.11) represents only one 
integral equation for e^ , since e~ is simply related to e"*" , and the equation for e"'^ is not an integral 
equation, as e""^ does not appear on the right hand side. 
At fixed A and T = 0, we get from ( p.ll| ) 

which are the order one excitations found in [^ for the ferromagnetic chain. The free energy per 



site in this limit is / = eo, which is the ground state energy found in 

The linear term in the specific heat C is obtained from the low temperature expansion of the 
free energy, or the entropy (3.5): 

C^-rg^Tg. (3.13) 

We find the linear term in the specific heat by computing the 0{T) term in the low temperature 



entropy [30, 31, ^, 15]. As T ^ 0, e^{X)/T scales as l/(Tcosh6A), which vanishes and gives no 
contribution to the integral (|3.9|) except when A ~ O(glnT). We rescale the equations ( |3.11 ) by 
making the change of variables A ^ A — | In T and consider the equations ( |3.11 ) at large A and 
small T. We define, at this range of variables, 4>^{X) = e^(A — | In T)/T. The integral equations 
become 



r(A) 



12e^ 



-6A 



iCi * In (1 + e""^ )(l + e-'^ ) - iCa * ln(l + e 



Ki * In 



(l + e-<^^)(l + e- 



K2 * ln(l + e 



,^""(A) ~ K2*ln (l + e--^ )(l + e 



(3.14) 



Differentiating ( |3.14 ) with respect to A, 

-2 X See"*^^ + Ki* 



dX 

d(j)- 
~dX 

~dA~ 



A'+ 



1 + ef 



+ 



l + e^ 



+ K2* 



l + e-^ 



Ki * 



+ 



-K2 * 



1 + e<^"^ 1 + e"?^ 



+ i^2* 



l + e? 



1 + 6-?^+ 1 + 



(3.15) 



where (j)'^ = d(p^ /dX. 

We rescale equations (3.1) in the same way. Let Pt{X) = Pt{X — ilnT), and recah that 



./3-n^ 



vr 



pf/(l + e^'^/^). Then 
ft+(A) = 
PtW = Ki* 

PTW = K2* 



^2V- + i^,*( ^* 



Pt 



+ 



l + e<? 
Pt" 



+ 



Pt 



1 + e<^+ l + e4 
Pt , Pt" 



l + e<^ 



Ko * 



Ko * 



/5r 



l + e? 



~ns 



1 + e-^"' 



+ 



■ 1 + e<^+ 1 + e<^ 
Comparing equation ( |3.16| ) to ( 3.15| ), we see that 



Pp 



k± 



1 



T 



1 



&TT dX 1 + e"^ 



67r dA 1 + e? 



(3.16) 



(3.17) 



The entropy can be evaluated in this hmit. The A — > 00 and A — > —00 hmits make the same 
contribution to S (and /). Therefore we write 



^ - 2 E /"dA{/5(^ln(l + e<^') + pfln(l + e-^')} 



2 X 



where 



/3=+,— ,ns 

7^ r /■</'''' (°°) 
67r L 7<^+(_oo) 

5(</') = 



</■ (00) 
(^-(-00) 



ln(l + e*) ln(l + e~ 



+ 



(/)"'' {-00) 



(3.18) 



(3.19) 



l + e'P ' l + e- 

The hmits (/)(iboo) are found from equations (3.14). In these limits, the integrals can be performed 
by taking the log out from under the integral sign, and integrating only the kernel. Let cj) denote 



10 



the asymptotic value under consideration (j){X 
equations 



±00). Then, for A = 00, we obtain a system of 



1 



1 



ln[(l + e'^ )(1 + e""^ )] - - hr(l + e 



1 



ln[(l + e-<^^)(l + e- 



Therefore the upper hmits of equation ( |3.18| ) are 

(/)+(oo) =(/)~(oo) =-ln2, 
At A = —00, (j)^ = 00, and 



--ln(l + e- 
^ln(l + e-^ 



(/.""(oo)=ln3. 



-ln(l + e- 



(3.20) 



(3.21) 



and thus the lower limits in ( p.l8| ) are 

0^( — 00) = 00, 0^( — 00) = 



In 



ri + \/5 



'(-cx)) =ln 



1 + ^5- 



(3.22) 



(3.23) 



We show in the appendix that how to express the integral ( |3.18| ) in terms of dilogarithms 
Using the identities described in the appendix, we find that in the low temperature limit, 

AttT 



S 



45 



and from (|l.5| ) we see that, with vp 

the well known conformal limit of the 3-state Potts chain (IT 



(3.24) 
This is the central charge of 



3 ||5| , the central charge c = | . 

33|, Isll computed in the limit (^) . 



This verifies that the limits ( |1.1D and ( |1.4[ ) smoothly connect together, and there are no additional 
length scales in the problem. 

In the calculation above, we find that, although in the zero temperature limit the ener- 
gies e~,e"'^ vanish, they contribute to the low temperature specific heat, i.e. the functions 
(/>~(A), i?i>"'^(A) do not vanish. This is a manifestation of the feature seen in [^, that although 
the energy can be expressed, as in (|3.8|), in terms of only p^, the number of states e with energy 



E depends on pf, and therefore on p~ and /O^^, as it depends on m-^mns in [^. In our case, 
the densities p~ , p^^ enter the expression for the free energy via the entropy S, which counts 
the states. Computation of the specific heat at low temperature depends sensitively on correct 
counting of states, as was seen in Q for the anti-ferromagnetic case. Note that from the counting 
rules of y], and from the equations (^.ij), the number of states near P = increases as the number 
of excitations increase. This causes the central charge to be larger than 1/2, the value we would 
expect if the excitations Pp were fermions. 



11 



4 Ant i- ferromagnetic chain 



To find the low temperature behavior of the anti-ferromagnetic chain, we rewrite ( ^.18 ) in term of 



the zb2s, ns densities, which we know to be the order one excitations for this hamiltonian |^, y, |3[- 



The equations ( 2.18 ) become 



±2s/\^ _ 3 ^ 17 , ( „^s , -2s\ -JJ , ns 



vr V2 cosh 3A =F 1 t^ ^ ^ 

_3 ^ 

TT cosh 3A 

where the kernels are: 



pTW = ^_^,, - K^*{pI' + pf') -2K,* p-^^, (4.1) 



-TT / , X 3 -TT / , X 6 cosh 3A , , , 

KiiX) = -, Ks A = -^ . (4.2 

'^ ' 27rcosh3A' '^ ' V2^cosh6A ^ ' 

Note that, from these equations, in the antiferromagnetic case the density of available states 
always diminishes with increasing particle densities. This was seen in [P, ^, where the number of 
available states decreases as a function of m2s,m_2s,rn,„s. 

The entropy in terms of these densities looks the same as (|3.5| ), but now we sum over the 
three independent densities a = ib2s,ns. We also express E in terms of the ib2s,ns densities. 
Minimizing the quantity ( |3.4D with respect to the three particle densities Pp"^ , p~'^^ , p'"'^ now gives 
the free energy in terms of these densities: 

/ ^ e-„ . r /°° .A{^'"<A+ '-""'"'^ + ^""' + --"•'^"^' + ""'^ + '-Pr^ }. (4.3) 
J-oc ^ 7r(A/2cosh3A-l) 7r(^/2cosh3A + 1) 7rcosh3A J 

where 

eo = eo + ^ / dX -^ = 3 ^ + - = 2.097 . . . , (4.4) 

vr2y_oo cosh 6A(^/2 cosh 3A-1) SVS vr 

and e" satisfy the integral equations 

e^'-W = ^^ ^3^ ^ ^ + t{Ki * ln[(l + e-"/^)(l + e-'"/^)] + K, * ln(l + e-"VT)}^ 

e-(A) = -^ + T{K2*lnUl + e-^''/^){l + e-^-''/^)] + 2K^.Hl + e-^"'/'')}. (4.5) 
cosh3A '- L J 

At T = 0, we see from (14.51) that 



,±2s _ ^ ,n. 



\/2cosh3A=Fl' cosh3A 



^r = ^r^, (4.6) 
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and from ( |4.3D the free energy is / = eo- These are the order one excitations and ground state 



energy for the anti- ferromagnetic hamiltonian found in |25, 



Again we compute the low temperature hmit of the entropy, but now we rescale the integral 
equations by changing variables to A ^ A — g InT. Defining </>^(A) = e(A — g lnr)/r we have 



/2s 



/-2s 



2 X ^e-3^ + t[Ki * ln[(l + e-^'')(l + e"-^ '')] + K2 * ln(l + e-"^"')} 

)(1 + e--^"'^)] + 2Ki * ln(l + e-*"')},(4.7) 



2 X 3e"^^ + tIk2 * In 



[1 + e-^'' 



so that 2s and —2s are symmetric in this limit, which was a feature seen in [^. Again, differentiat- 
ing ( |4.7| ) with respect to A and comparing to the density equations { \i.l\ ) rescaled as A ^ A— | InT, 
we see that 



4 



1 



P = 2s, — 2s,ns. 



Sir dX 1 + e't" 
The entropy is calculated as in ( |3.18| ). From ( |4.7| ) we find the limits (j)°'{±oo) 

(/,2^(-0o) = 0-2«(-CX))=0""(-Oo) =00, 



-2s, 



00 



</,2*(oo) =ln2, 



00 



In 3. 



The entropy is 



(4.8) 

(4.9) 
(4.10) 

(4.11) 

This, with vp = 3/2 [^], gives a central charge c = 1, which is the central charge of the conformal 
limit of the model pO| , p^ j. Again, this verifies that the limits (1.1) and ( |1.4| ) commute, and there 
are no additional length scales. We see that the fact that pf tends to decrease with increasing Pp 
causes the central charge to be smaller than the value 3/2 one would expect for pure fermions. 



Svr 



In 2 

2 / g 

00 



■In 3 

+ / 9 

00 



which, using the dilogarithmic identities in the appendix, gives 

27tT 



5 Q = ±l 

In the thermodynamic limit, quantities which are not of order M are irrelevant to the calculation. 
Therefore, we do not expect the value of Q to affect the thermodynamic equations. This is indeed 
the case. 

In 1^], the counting rules for Q = ±1 were found. It was shown that the counting rules 
depended on the value of the numbers m+-|- and m ^, where 



m+H 



m_ 



0,±1. 



(5.12) 
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Again we note that this difference is not of order M and we do not expect it to change the 
thermodynamic equations. The sum rules ( 2.16| ) for this sector are changed to: 



M-1 = m2s + 2m„s + 3m_2s +rn_+ +m++. (5.13) 

This is only different from the (5 = sector by a term of order 1, due to equation (|5.12| ) and the 
fact that m y = m_. Therefore to order M, the sum rules are identical to ( |2.16| ). 

For m+-|- — m y. = ±1, the completeness rules for Q = 1 are the same as for Q = Q. For 

m++— m I- = 0, there is spectrum doubling: The integers I^ are shifted from those of + by ±1/2, 

and those of —2s are shifted from /•" in the same way, both signs giving the same energy level. 
The shift does not affect the thermodynamic limit of equations ( |2.15| ). The spectrum doubling 
gives rise to an additive term of order one in the entropy (which counts the number of states). 
Since the entropy is of order M, again this term is not relevant in the thermodynamic limit. We 
conclude therefore that this sector is identical to Q = 0. 

In [^, the difference in counting rules for the sectors Q = ^1 gave rise to different branching 
functions in the modular invariant partition function from those of Q = 0. However, each term 
in the modular invariant partition function gives the same specific heat. Due to the modular 
invariance property. The specific heat is found from the limit (7 ^ 1 of the partition function 
of IP]. However, the partition function is invariant under modular transformations, where, if 
q = exp(27rzr), the transformation r — > — l/r leaves the partition function invariant, so the specific 
heat is obtained from the g — > limit of the partition function. The same transformation sends 
each branching function into a linear combination of all other branching functions. Therefore, 
each branching function has the same g — > 1 behavior. In the thermodynamic calculation we do 
not see this difference between the sectors (5 = and (5 = ±1. 

6 Discussion 

In H, the order one excitations ( [4. 61 ) were used to compute the partition function of the anti- 
ferromagnetic chain in the limit (I]^) . This is identical to the modular invariant partition function 
of the conformal limit of the model ||l^, ^ and so gives the same specific heat as ( [4.11 ). For 



the computation in [^, it was necessary to carefully consider the information about the way the 
maximum integer I^ax of the finite size system in ( ^.15 ) changed with the number of particles 



in the system. This is because those integers correspond to the energies close to zero, e(P) ~ 0, 
or in the language used here, the large A behavior of the energies, which is the region which 
contributes to the specific heat. The maximum integers decrease, at finite M, as more particles 
are added to the system. In Q this phenomenon was referred to as an infrared anomaly. This 
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represents counting rules for the excitations beyond the fermionic exclusion rule, and is the phe- 
nomenon responsible for the central charge being different from that of fermions. In the case 
of the anti- ferromagnetic chain discussed in [0], the infrared anomaly was repulsive: fewer states 
were available as the number of excitations was increased than would be available for fermions. In 
the case of the ferromagnetic chain, both repulsive and attractive infrared anomalies are present, 
but the total infrared anomaly is attractive. 

In the thermodynamic limit, we discard the information about the maximum integers /^. 
Nevertheless, the integral equations (f4.5D contain the information about the way the density of 



available states depends on the particle density, represented by the density equations (4J). This 
information enables correct counting of states, using the entropy. We see that this information 
gives the same specific heat as the counting of ref Q, but, as we saw in section 5, does not show 
the difference between the different Q sectors. 

When computing the low temperature specific heat, the region of A which contributes to the 
free energy as T — > is the A ~ |lii2^ limit in the ferromagnetic case, and the A ~ |lnT 
limit in the anti-ferromagnetic case. This limit is the P ^ limit [^, which corresponds to the 
lowest lying order one excitations in |^] . The careful counting of states there is paralleled in the 
computation here by the rescaling of the integral equations and the density equations in the limit 
T — s- and A — > oo. 

Finally, we note that these computations are related to the thermodynamic Bethe ansatz 
method of |Q , a point which is discussed in some detail in 1 22 1 . 
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A Appendix: Expression of entropy via dilogarithms 

The Rogers dilogarithm is defined as 



^^■^^-\i:<'^^'my 
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Making a change of variables in the expression for the entropy ( |3.18| ) to / = 1/1 + e*^, the 
entropy ( |3.18| ) is expressed in terms of L{x) [p^]: 

So.-^ (2L{^) + L{\) - L(l) - 2L(^^)] . (A.2) 



We use the identity on Rogers dilogarithms |l32| : 



g / s.nV/„) \ 2(n-j) 

where L{1) = 7r^/6 |3g]. Using this identity with, n = 6, we find that 

2L(l/3)+L(l/4)=L(l), (A.4) 

and, with n = 5, we see that [pq]: 



'3-^' 



.vr 



2 



'M^^ rv (^■^' 



Therefore equation (A_^) gives 



In the anti-ferromagnetic case, we again use ( |A.4| ) with n = 6, and entropy is: 

AT AT IttT 

S ^ ^(2L(l/3) + L(l/4)) = — L(l) = ^-. (A.7) 
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